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Water drops are formed at cloud base by condensation upon nuclei and as they rise they 
grow by vapour diffusion in the slightly supersaturated environment of the clouds. Turbu- 
lent mixing between cloudy air and undersatured air entrained from outside produces fluctu- 
ations in supersaturation, not linked to changes in vertical velocity, which cause broaden- 
ing of the condensate spectrum and the rapid production of droplets large enough to engage 
in growth by coalescence. The probabilities of permanent union or the production of satel- 
lite droplets following the collision of a pair of raindrops is a sensitive function of 
several parameters. In some circumstances electrohydrodynamic bursting may influence the 
the properties of clouds. 


1 -Introduction 

Liquid drops produced by natural processes in the Earth’s atmosphere vary in size from 
around 0.1 urn to 1cm — a range of about twelve orders of magnitude in volume. Thus it is 
not surprising that such drops are involved in a large number of physical processes C'ld 
phenomera; some of which are of crucial importance to large-scale energy transport, the 
cleanliness of the atmosphere and the fertility of cultivated land. It would be beyond the 
scope of this short article to attempt a detailed review of each of these processes. Instead 
they are treated, in* Section 3, in a cursory manner which, it is Hoped, will permit their 
range and their interdependence to be outlined} and one process - the diffusive mixing of 
dry and cloudy air, which is considered to be fundamental to the development of natural 
clouds - is selected for more comprehensive discussion , in the following section. This 
particular topic has been chosen in view of recent concentration of attention upon it by 
several groups of scientists. 

2 . The diffusive mixing of dry air with a population of cloud droplets 

The sub-adiabatic liquid-water-contents generally observed within clouds from which ice 
is absent - or present only in insignificant quantities - are a consequence of the dilution 
of the clouds, during their growth, with environmental, undersaturated air. Recent work 
{Baker and Latham, 1979^ Baker, Corbin and Latham, 190C^ Telford and Chai, 19B04 Telford 
and Wagner, 1991^ has indicated that the non-uniform manner with which the environmental 
air mixes with the cloud gives rise to fluctuations in supersaturation - not linked to 
changes in updraught speed - which cause spectral broadening. A small fraction of the 
droplets move through the condensate spectrum to the coalescence ’ge of growth several 
times faster than predicted by classical descriptions of the ent *^ent process (for 
example, Warner (1973 )i Mason and Jonas, (1974)^ Lee and Pruppach ^977)'), in which the 
mixing between environmental and cloudy air at any level was assumti to be instantaneous 
and uniform. The droplet spectra predicted by Baker et al agreed well with those measured 
in cumulus by Warner (1969). 

Baker et al did not specify the nature of the entrainment process but assumed that a 
volume of environmental air mixing instantaneously with cloudy air at a given level reduced 
the droplet population there by means of a combination of dilution and total evaporation 
of a fraction of droplets of all sizes. Telford and Chai assumed that outside air entrained 
into the cloud at its top became negatively buoyant as a result of droplet evaporation, 
and then descended, as a saturated stream, to mix with and thereby dilute cloud air at a 
lower level. The common consequence of both these mixing processes is a droplet spectrum 
identical in shape to that prior to mixing, but with reduced number concentrations in all 
size categories. Subsequent ascent of the mixed volume produces an enhancement in super- 
saturation as a consequence of reduced competition for water vapour; as a result, the drop- 
lets grow more rapidly than those in adiabatic regions of the cloud. 

In neither of these studies was account taken of the finite rates at which cloudy and 
environmental air (or adjacent regions of clcud, with different properties) will mix. In 
this paper we present a model of the turbulent mixing of a spherical blob of droplet-free 
air, of original diameter Xg, temperature Tg and supersaturation S (< 0), and a cloud, 
represented by a spherical shell surrounding th^ lob of outer diameter Xq, temperature Tq, 
supersaturation 5 (»0) and liquid water content L. The blob contains a distribution nCmg) 
of-cloud condensation nuclei, comprised of NaCl particles of mass m 3 in equilibrium at the 
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relative humidity H « S 4 I . This distribution is identical to that aasumsd to be contained 
within the droplets in the surrounding cloud, so that when turbulent mixing occurs at any 
level the total concentration of particles (droplets plus CCN) remains constant. The entire 
system is assumed to be moving upwards with a speed W. 

The turbulent mixing of all properties transported is represented by diffusion with a 
single coefficient K determined on dimensional grounds from and e, the rate of dissipa- 
tion of turbulent energy within the cloud; K . If we assume that the turbulent 

motions transport energy, liquid-water, water-vapour and dry air simultaneously, we have 


|t(R,t) ^ ^^(R.t) _ K7^h(R,t) - qf,(R,t). 

(1 ) 

ffl(fi.t) ^ u^(R.t) - Kv^o,(R.t) - Q,(R,t). 
at ^ ^ 

(2) 

l^(R.t) ^ y^(R,t) _ _ Q^(R,t) 

(3) 

^(R.t) ^ y^(R.t) _ (<„2o^(R,t) - Q^(R.t) 

(4) 

^(r^.R.t) ^ y|Q(r, .R.t) _ ^p^t) « q (r.,R,t) . 

ot pz 1 * n 1 

(5) 


P , and D are the densities of dry air, liquid water and water vapour, respectively; 

Eft) the number density of droplets with radii in the interval ( r , ri+dr^ ) ; and 
£ is the redial vector. The energy term 




( 6 ) 


where and C, are the specific heat of dry air and liquid water respectively and 
Ly(»2.2oKJ g**'; is an average value of the latent heat of vapourization of water. Equation 
( 5 ) defines the temperature T in terms of the parameters of the problem. 


The source terms Q are; 

Q^(Eft) « -rt^Q^(Ef t)+Pg(R, t)gW(t) 3 
Q„(fift) • -4np lii(r . ,R, t)r?f . (£, t) 
Q, (R.t) - -Q„(R,t) 

^ ''n(ri.R.t) 

- P ( R.t) • 


Q^(R,t) - -P^( 


R t) ami) 

|T(R,t) ^ R^T(R,t 


(T) 

( 8 ) 

(9) 

( 10 ) 

( 11 ) 


where is the density of liquid water, g the acceleration due to gravity and the gas 
constant per kilogramme of air. 


Qh(E.t) 

' CC^P^CR.t) + Cj(Pj^(R,t) + P^(R,t))J 


The droplet growth equation is 
1 

where 


DP (T)/P r . Bm 

■ (r.-fs) S(a,t) - 7 3 

X 1 


( 12 ) 


(13) 


(14) 


P,(R,t) « 1.B X 10^exp(-5400/T(R,t)) 

is the saturation vapour density (g m*^), D is the coefficient of molecular diffus^wo 
(cm^s*M, the supersaturation S is expressed as a percentage, r. is measured in micrometres, 
mg in grammes, and the constants A and B possess the values 0.1 i5 and 1.4 x 10^^ respectively. 
aT*5um) is a characteristic length associated with the non-ideality in water condensation. 
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The rate of change of auparsaturation 


4(R,t) - (100 + s(R,t)) 






5,4 X lO'^Q.^(R.t) 
T^(R.t) 


3.42 X 


T(fl.t) 


10~^W(O 


(15) 


The precieion of the foregoing equatione could be increeeod» but such refinements were 
deemed unnecessary, in view of the very approximate relationship between this diffusion 
formulation and the mixing process within clouds. We note that the diffusivity representa- 
tion is best applied to situations in which exceeds the characteristic spatial scale of 
the turbulent eddies. Although buoyancy gradients generated by evaporation will deform the 
volume under examination, we assume that spherical symmetry is preserved during the mixing 
process . 

Figure (1) illustrates the predicted time variation of the fields of liquid-water-content, 
if supersaturation, 5, and temperature, T, resulting from the mixing of a blob of air of sire 
Ig • 100m and relative humidity 80% with a cloud of sire X • 200m and liquid water content 
L • 0.5^ m"3, W ■ 0 m s" . In the first case illustrated (S,B,C) the cloud consists of 
droplets of radius 3um and in the second {D,E,F), r • 20um. Since the turbulent energy diesis 
pation rate e • 10"^m^s~3 the characteristic mixing times 


Yj- » (X^^/e)^ ■ 1 60s 

(16) 

TB " (Xg^/e)^ • 100s 

(17) 


We see, when r * 3um, that evaporation of droplets carried by turbulence into the interior 
of the blob is a very rapid process - presumably because of the high surface/volume ratio 
of the droplets - so that the supersaturation rises from -20% to close on 0% in about 10s, 

The temperature difference (originally 2K) follows a similar pattern - since its changes are 
governed by evaporation when the latter is rapid - lagging just slightly behind changes 

in 5. Evaporation is confined more-or-less entirely to the region originally occupied by 
the blob. The liquid-water-content distribution takes considerably longer to level out, 
because evaporation inhibits its increase in the early stages of mixing, but nevertheless 
the gradient has been substantially reduced after 20 seconds, a period much leas than the 
turbulent mixing times and - 

rigure 1 shows that when the cloud consists of larger drops (r ■ 20um) the rate of 
increase of minimum supersaturation 5 towards zero is reduced - because evaporation j.a less 
effective - but is still essentially completed within 20 seconds. The temperature equaliza- 
tion curves again follow those of S, and there is a corresponding increase in the rapidity 
with which the gradient in L levels out. In this case turbulence distributed the under- 
^eturation significantly into the cloud. Mixing is largely completed after 20s. 

Figure (2) shows the time-development of the initially monodisperse spectrum (rw20umi the 
droplet concentration IM ■ 1 5 cm*^ ; L » 0.5g m~^ ) at R • location of the original 

interface between the cloudy and undersaturated air (H ■ 00%), W m 0.0m s*^ . It is seen 
that the combination of evaporation and turbulent mixing introduces many new categories into 
the droplet size distribution. 

Figure (3) presents the calculated spatial fields of L and S, at various times, for two 
different values of e, the turbulent energy dissipation rate, W ■ 0 m s* . In this case^ 
the original spectrum, illustrated (A) in Figure 5, has a liquid water content L • I.Og m" . 
It contains ten size classes (radii between 3 and 12.5am) and is ba»d on ones encountered* 
in field studies on thu summit of Great Dun Fell (Blyth et al (190On, In this case Xjc ~ 
130s when e ■ 1 CJ®m^s“** ( A , B ) and XtC ^ 270s when e ■ 1 Cl"^m^s*^ ( C , D ) ; the corresponding valuea 
of Xj 3 are 74s and 160s. When the mixing is rapid (A,B) the supersaturation becomes uni- 
formly distributed within 20s, and gradients in L are eliminated well before 50s, With 
slower mixing the equalization of both L and 5 are slower, but are completed in times less 
than tjq. In both cases undersaturations are created within the region occupied by the ori- 
ginal cloud; a characteristic feature of the situations in which the cloud spectrum contains 
large drops - with correspondingly high values of evaporation time (Baker et al,(190Or). 

In figure 4 we present, for the cases covered by Figures 1 and 3, information on the rates 
at which the total amount of evaporation and the distribution of liquid water approach their 
final values. f(in curves A,B,C,D) is the ratio of the total evaporation that has occurred 
at time t (expressed, in Figure 4 as the dimensionless time t/T^) to that which has occurred 
when the cloud and blob are ^ully mixed. In curves E,F, G,H 
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A 


o't.O 


1 - AL/L 


(18) 


where AL • (Lq - Lo)+ is the difference, at time t, between the liquid water contents in 
the centre of the original blob (R«0) and at the edge of the cloud {Hm\^/2)i L is the ori- 
ginal liquid*water-content in the cloud. 


The curves displayed in Figure 4 show clearly that: evaporation is completed long before 
the gradients in L are eliminated; uniformity in L is achieved - generally - in times much 
less than evaporation proceeds most effectively when droplets are small and mixing is 

fast. Further calculations (not illustrated) show that uniformity in S is achieved more 
quickly for increased values of L. 


In Figure 5 we present various spectra resulting from the mixing with undersaturated air 
of a volume of cloud possessing an original liquid-water-content L ■ 0.5g m“^ and the sire 
distribution (A) mentioned earlier, W *0.0 m s"*^ . For the particular dimensions chosen 
(Xc « 144m, Xb • 60m) the final - well-mixed - value of L, if the humidity of the entrained 
air is 70?^ and the temperatures are as indicated, is Q.32g m*^. Curve B presents the 
classical spectrum, with L ■ 0.32g m*^, which was obtained by the application of the drop- 
let growth equation to all droplets in the original spectrum. All droplets have evaporated, 
so the peak radius has shifted substantially (whereas the peak concentration is virtually 
unchanged) and the reduction in radius is greater for the smaller droplets. Curve C, L ■ 
D.32g m*3, is obtained from our diffusive model, with c ■ 10“3m2g-3^ This spectrum is seen 
to have a peak-radius only slightly below that of the original spectrum (A), but a substan- 
tially reduced peak-concentration. Thus the spectrum C, predicted from our model, lies 
between the classical spectrum B and the extreme inhomogeneous spectrum, but is very much 
closer to the latter - which is not shown since its similarity to B would cause confusion 
in presentation. 

In Figure 6 we present three spectra resulting from the turbulent mixing of spectrum A - 
the same as in Figure 5 - with a blob of undersatureted air, W • 0 m s* ' . The dimensions 
and temperatures, given in the legend, are the same for all three calculations. In B and 
C the relative humidity H ■ 70% and the final water content L • 0.32 g m*3 . In B, 
e • 10“^m^s“^ and in C,e ■ lO’^m^s*^. Curves B and C are very similar, but B is more 
’inhomogeneous* - smaller reduction in peak radius. - consistent with the associated slower 
mixing in this case. Curve D (L * 0.1 8g m*^) was obtained with H » 40% and € • l0"3m28-3. 
This spectrum, although very different from one predicted by the classical description of 
mixing/evaporation (see, for example, Figure 4, Curve B), is substantially different from 
the equivalent ’extreme inhomogeneous’ or ’dilution only’ spectrum (again not showni, in 
which th3 shape of curve A would be preserved, and the mean radius would be the same as in A. 
This large departure from the extreme inhomogeneous spectrum may be attributed to the greater 
role of evaporation (relative to dilution) in producing the final spectrum (D). 


In Figure 7 we show the effect of doubling the radius of all droplets in the spectrum A 
illustrated in Figure 5 and 6, whilst maintaining the liquid water content L » 0.5g m*3^ 

W ■ 0 g m*3, jhe final spectrum (C) after evaporation and mixing is closer to the classi- 
cal one B than in the earlier cases, presumably because for these large droplets the ratio 
of evaporation to mixing times, (Baker et al (1980)2) is no longer much less than unity. 

In Figure B we attempt to examine, for a wide range of original conditions, the extent 
to which spectral changes resulting from turbulent mixing of undersaturated and cloudy air 
conform to the classical (homogeneous) picture or to the alternative extreme discussed by 
Telford and Chai, and by Baker et al. For ease of interpretation, these calculations were 
made with originally monodisperse spectra. 


On the classical picture, to a first approximation, the concentration of droplets, N, 
in an originally monodisperse spectrum, is unchanged from the original value N^, as entrain- 
ment reduces the liquid water content from L- to L. Or. the extreme inhomogeneous picture, 

N is reduced in proportion to the reduction in L, (ie N ■ (L/Lq)Nq). On the classical 
picture the water contert is reduced in proportion to the reduction in the radius r(ie 
r3 « (L/Lo)ro» where Tq is the original radius). On the extrame inhomogeneous picture r 
is unchanged as L decreases. Thus we may define a parameter 


9 m (N/N )(L/L )/(R/R„)^ 
0 0 o 


(19) 


which possesses a value qp « 1 for all L/L^ on the classical description of mixing and is 
given by qp • (L/Lq) 2 in the extreme inhomogeneous case. 

These two extreme curves are plotted in Figure 6, togethex with individual values of qp 
emanating from a variety of calculations, using the diffusive mixing model. 
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The most important reault la thatf for all situationa conaidarad, tha axtrama Inhomo* 
gan eoua daacription of turbulent entrainment pro v idea a much batter approximation 

to reality than doaa tha claaaical daacription • Tha departuraa from can 

be quite aubatantial, however. We see from Figure 6 that this departure incraaaaa withi 
decreasing increasing H (at constant L/Lo){ dacreaaing L| decreasing cf and (proportion- 
ately), decraafing L/Lg. All these tendencies are consistent with tha basic idaaa» advanced 
by Baker at al| that { T ) the time constant ratio which equals zero in the axtrama 

inhomogeneous case, increases with decreasing values of L and X and with increasing c 
and and (2) that as (L/Lq) is reduced (by reducing H) the contribution of dilution 
to the overall reduction in water content is reduced. 

3. QihSL .MggMttgj 

It is inappropriate, in an article of this limited length, to do other than to list further 
important processes, occuring in clouds* which involve drops. For a comprehensive and author!* 
tatiye treatment of these the reader is referred to the excellent book of Pruppacher and Klett 
(t97&rVThe processes include: condensation via heterogeneous nucleation; collision; coales- 
cence; satellite drop production; electrohydrodynamic bursting; mechanical disruption; 
scavenging of particulate and gaseous pollutants; supercooling; freezing; accretion; 
secondary ice particle production; splashing; and lightning initiation by cc^^ona. 
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Figure 3: Calculated radial distributions of liquid-water-content , L, and superaaturation 

S at various times t during turbulent mixing of a water cloud with a, spherical 
blob of undersaturated air centrally embedded within it at t«0. Ten cotegory 
spectrum. 60m; ^ «144mj| m. « L ■ I.Og m”^; H ■ 1D%*, Tu« 261k; 

T^- 279K. A, B, e « C, D, c - IQ-^m^s-^;. 

1, t « D.Os; 2, t« 8s: 3, t « 20s: 4, t ■ 50s. 
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Figure 5: Droplet size drstribution (C) calculated by application of the diffusive mode^ 

to the mixing of a apherical blob of undereaturated air with a water cloud pos* 
seasing the size distribution fA) at time t » 0. Spectrum B is calculated 
on the classical model of mixing. The liquid water content in A is 0.5g m*^ 
and in B and C 0.32g m"^. Xo ■ fiOmi X- • 144m{ €• 10*"^m^s“^l H ■ 70lt| 

T0-2B1KJ Tc« 279KJ m, • 10-^5g. 
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i Droplet tare dietribution (C) calculeted by application of the diffutiva modal 
to the mixing of a apherical blob of undaraaturatd air with a water cloud pot- 
tetaing the tire diatribution (A) at time t • 0. Spectrum (B) it calculated on 
the clcatical model of mixing. The liquid water content in A ia 0.5g m"* and 
and in B and Z pi3?g m"3, Xp* 60 m; Xc* H4m> c« 1p-2m2#-3| h ■ 70< Tp«2B1Ki 
Tj.-279Kj m,-10-^^g. ® 
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